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Using the kinetic theory we discuss how the particle and energy densities 
of the quark-gluon plasma fluctuate in a space-time cell. The fluctuations in 
the equilibrium plasma and in that one from the early stage of ultrarelativistic 
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the parameters involved the fluctuations appear sizeable in both cases. 
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The quark-gluon plasma, which is expected to be produced in energetic heavy-ion 
collisions, is usually described in terms of the averaged quantities like baryon or energy 
density, temperature, etc. However, it is often important to know how these quantities 
treated as local ones fluctuate around their average values. For example, it has been 
argued in our papers [IJ that the color current fluctuating around zero can initiate the 
plasma instabilities. Gyulassy, Rischke and Zhang have recently shown that the 
hydrodynamical evolution of the quark-gluon plasma can be significantly influenced 
by the energy density fluctuations at the initial state. 

In this note we derive, using the kinetic theory methods, simple analytic expres- 
sions describing how the particle and energy densities fluctuate in a space-time cell 
of the volume Ax 3 At. Then, we use the derived formulas to estimate the particle 
and energy density fluctuations in the equilibrium plasma and in the parton system 
from the early stage of ultrarelativistic heavy-ion collisions at RHIC or LHC Since 
the approach, which is employed, is rather simplistic our results should be mainly 
treated as a guiding for a more quantitative study in the future. 

Let us start with the fluctuations of particle density p(x), where x = (t, x) is a 
four-position. The average density is expressed through the distribution function / 
as 

MO) = (P) = I /(P) > W 

where p is the particle momentum. The system is assumed be on average homoge- 
neous and stationary and consequently, the distribution function is independent of x. 
The momentum distribution is arbitrary. 

The density correlation function for the classical system of noninteracting particles 
is given by a well known formula, see e.g. 0, as 

A(x) ^ (pixM^)) ~ (P) 2 = / /(P) ^ (3) (x- V) , (2) 

where x = X\ — x 2 = (ti — t 2 , xi — x 2 ) and v p is the particle velocity equal p/E p with 
E p = i/p 2 + m 2 - Due to the average space-time homogeneity the correlation function 
A(x) depends on the difference of Xi and x 2 only. The space-time points (ti,xx) and 
(t 2 , x 2 ) are correlated in the system of noninteracting particles if the particles fly from 
(£i,Xi) to (t 2 ,x 2 ). For this reason the delta 5^ (x — vt) is present in the formula (Q). 
The momentum integral of the distribution function simply represents the summation 
over particles. Applicability of our classical approach to a quantum system such as a 
quark-gluon plasma is considered below. We also discuss there when the partons can 
be taken as noninteracting as in eq. (||). 

We get the particle density fluctuations in a given space-time cell averaging the 
correlation function A(x) over the cell volume. Specifically, 

(p 2 )-(p) 2 = Jd 4 xA(x)A(x), (3) 
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where the coarse-graining function A(x) is chosen to be of the gaussian shape i.e. 

1 / t 2 x 2 

A(x) = , 9 t- t , exp( 



47r 2 At Ax 3 *V 2 At 2 2 Ax 2 

with At and Ax denoting the root mean square of the cell size in time and space 
direction, respectively. 

Substituting the correlation function (Q) into eq. (|3]) we get 

In 2 ) - In) 2 = I / IM ( 4 ) 



As seen, the fluctuations remain finite for At = but grow to infinity when Ax — > 0. 

For massless particles v 2 = 1 and the formula (|j) essentially simplifies. The 
relative fluctuations then are 

(p) 2 (2tt) 3 / 2 Ax 2 VAx 2 + At 2 (p> ' 1 ; 

and depend only on the average density and the size of the space-time cell. In the 
case of massive particles, a similar simplification occurs for At = 0. 

One easily generalizes the above considerations to the energy density fluctuations. 
The average energy density and the respective correlation function read 

< £ > = /(&^/(P). (•> 

W(x) S (e( Xi )e(x2)} - (e? = J El /(p) S<V(x - v,t) . 
One finds the energy density fluctuations in the space-time cell as 

d 3 p E 2 /(p) 



(2n) 3 / 2 Ax 2 J (2vr)3 ^ Ax 2 + v 2 At 2 

For massless particles the relative energy density fluctuations are 

(s 2 ) ~ (s) 2 1 (e 2 ) 

(e) 2 (2tt) 3 / 2 Ax 2 ^ Ax 2 + At 2 (e) 2 ' 

with 



(7) 



One sees that in contrast to the relative particle density fluctuations given by eq. (|5|), 
the energy density fluctuations (|7|) depend not only on the average energy density 
but on the on the energy second moment, i.e. (e 2 ), as well. 
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As already mentioned, our approach is classical. There are two quantum effects: 
particle localization and bosonic or fermionic statistics which should be considered 
before the approach is applied to a quantum system. If the cell size Ax is much larger 
than the average length of particle de Broglie wave A, i.e. 

Ax > A , (9) 

the particles can be treated as well localized. One also finds analogous requirement 
for At, but for the ultrarelativistic system under consideration, the space and time 
scales, as the momentum and energy ones, are usually of the same order. Thus, 
we further discuss only the condition (|^), which is checked below for the two cases: 
the equilibrium quark-gluon plasma and the parton system from the early stage of 
ultrarelativistic heavy-ion collisions. 

One takes into account the effect of quantum statistics substituting the distribu- 
tion function / by f(l ± /), where + is for bosons and — for fermions, in eq. @. For 
the equilibrium quark-plasma the numerical effect of the statistics is small because 
the quark (fermionic) and gluon (bosonic) corrections, which are of comparable value, 
act in the opposite directions. In the case of the parton system from the early stage 
of ultrarelativistic heavy-ion collisions, where gluons dominate, we also neglect the 
effect of bosonic statistics. The correction is expected to be small but significantly 
complicates the computation. The point is that the parton density is large at the 
early stage of the collision but the phase-space density, which matters for the quan- 
tum statistics effects, is reduced due to the large longitudinal momentum range of 
the partons. 

The quarks and gluons are taken as noninteracting in our approach. More specifi- 
cally, the partons are assumed to follow the straight line trajectories in eq. (0). Since 
the QCD color forces are of long range, the quarks and gluons always interact in a 
many-parton system. However, the effect of interaction is minor at the space scale 
which is much smaller than the inverse momentum transfer due to the interaction. 
Therefore, our formulas @ and (0) are basically correct as long as 

Ax < - , (10) 
Q 

where q is the characteristic momentum transfers discussed below. 

To obtain our final formulas @ and (|7|) the partons have been assumed to be mass- 
less. This assumption is correct for the case of equilibrium plasma at the sufficiently 
large temperature when the interaction can be treated as a small perturbation. The 
situation is less clear when the early stage of nucleus-nucleus collision is considered. 
Then, the partons are off mass-shell. Nevertheless we still treat them as massless i.e. 
we assume that the parton mass ('offeshellness') is smaller than 1/Ax. We adopt 
this assumption keeping in mind that the partons with the 'offshellness' [i decay after 
jjT 1 . Thus, they can be treated as noninteracting at a much smaller scale only. 

After these comments let us discuss the equilibrium quark-gluon plasma. The 
distribution function then reads 
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where (3~ l = T is the temperature and g g = 16 and g q = 24 are the numbers of the 
internal degrees of freedom of quarks (of two flavours) and gluons within the SU(3) 
gauge group. The plasma is assumed to be baryonless and the quarks, as gluons, are 
massless. 

Substituting the function ( |TTD to eqs. (|IJ,^) and (Q) we get, respectively, 



(p) = 


34C(3) T 3 


= 4.14 T 3 


(e} = 


377T 2 rp4 

30 


S 12.2 T 4 


e 2 ) = 


462C(5) jn 5 

IT 2 


^ 48.5 T 5 



with ((z) being the zeta Riemann function; £(3) = 1.202 and £(5) = 1.037. The 
relative fluctuations of particle and energy densities (||) and (0) then are 



(P 2 )-(P) 2 „ 0.124 



\ (p) 2 Ax 3 / 2 T 3 / 2 



(e 2 )-( £ ) 2 0.144 . . 

\ (e) 2 Ax 3 / 2 T 3 / 2 ' 1 ; 

where we put for simplicity At = 0. 

We estimate the average length of particle de Broglie wave as the inverse average 



momentum, which in turn is defined as y((p 2 )), where ((...)) denotes averaging over 
particles. Therefore, 



A 



(p) „ 0.292 



and the condition (Qj gets the form Ax T>0.3. 

In the equilibrium quark-gluon plasma the characteristic momentum transfer can 
be identified with the Debye screening mass ran, which for two flavours and three 
colors equals gT with g being the QCD coupling constant, see e.g. [|]. Then, the 
condition flID"D is AxT <^l/g. 

One sees that the two conditions, which allows us to treat partons as free and 
classical, can be fulfilled simultaneously if g^ 1 ^> 1 i.e. when the plasma is weakly 
interacting. Unfortunately, at the temperatures of order of a few hundreds MeV g 
is not much smaller than unity. Therefore, Ax should be close to T~ l that the two 
conditions are not badly violated. 

The result (0) can be helpful in choosing the physically reasonable size of the 
elementary fluid cell in the numerical hydrodynamical calculations. On one hand, 
the cell should small enough to get details of the density profiles, but on the other 
hand, the cell size should be sufficiently large to reduce the fluctuations. One sees 
from eq. (0) that at T = 200 MeV and Ax = 1 fm the energy density deviates 
from the average value by about 14%. The fluctuations increase to 41% for Ax = 0.5 
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fm. However, such a small cell is at the border line of applicability of our classical 
approach. 

The size of the space-time cell is sometimes dictated by the characteristic scale 
of the phenomenon under considerations. When the deconfmement phase transition 
is discussed, the cell size Ax should be identified with the inverse confinement scale 
parameter Aqcd — 200 MeV. At the critical temperature T c = 150 MeV the en- 
ergy density fluctuations are about 22%. Such sizeable fluctuations can significantly 
speedup the process of hadronization. However, one should keep in mind that near the 
phase transition the quark-gluon plasma is no longer a system of weakly interacting 
particles and the formulas derived above can provide only a very rough estimate. 

Let us now consider the nonequilibrium plasma from the early stage of ultrarela- 
tivistic heavy-ion collisions. It is commonly believed that the perturbative processes, 
which are under theoretical control, play an important, if not dominant, role in these 
collisions, see e.g. ||. We adopt this conventional point of view and discuss density 
fluctuations of hard and semihard partons which are characterized by the transverse 
momentum being relatively large, at least 2 GeV, when compared to the QCD scale 
parameter Aqcd- 

We estimate the particle (mostly gluons) density taking the numbers from ||, 
where it has been found that about 570 perturbative gluons are generated at the early 
stage of the central Au-Au collision at RHIC (\/s = 200 GeV per N-N collision) and 
8100 at LHC (-\/i = 6 TeV per N-N collision). Assuming that all these gluons appear 
in the cylinder of the volume nr^A 2 ^!, where r$ = 1.1 fm, A = 197 and / = 1 fm, we 
get the average densities 



(P) 



4.4 fm" 3 for RHIC , 
63 fm" 3 for LHC , 



which immediately provide (via eq. (||)) the density fluctuations. 

To estimate the energy density fluctuations one needs the parton momentum dis- 
tribution. We take it in the form which corresponds to the flat rapidity distribution 
in the interval (-Y, Y) i.e. 

f( P ) = -Le(y - y)e(Y + y ) h{ Px ) -\- , (is) 

2Y p± chy 

where y and p± denote the parton rapidity and transverse momentum. We do not 
specify the transverse momentum distribution h(p±) because it is sufficient for our 
considerations to demand that the distribution ( |l~3l) is strongly elongated along the 
beam axis i.e. e Y ^> 1. 

The QCD-based computations, see e.g. 0, show that the rapidity distribution 
of partons produced at the early stage of heavy-ion collisions is essentially gaussian 
with the width of about one to two units. When the distribution ( |i~3"D is used to 
simulate the gaussian one, Y does not measure the size of the 'plateau' but rather 
the range over which the partons are spread. If one takes the gaussian distribution 
of the variance cr and the distribution (O) of the same variance, then Y = V3 a. 
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One computes the average energy @ and energy squared (§) densities with the 
distribution ( |T3| ) as 

(e 2 ) = Sh2 \y 2Y (P1)(P) = ^W) 2 (P), 



and gets the energy density fluctuations 



(g) - (e) 2 0.178 
^ (e) 2 Ax 3 / 2 




(14) 



where At = 0. 

Taking the values of F given in § i.e. F = 2.5 at RHIC and F S 5.0 at LHC 
and the gluon density, which has been estimated above, one gets the energy density 
fluctuations. Eq. ( |14"D tells us that the fluctuations at RHIC are 13% for Ax = 1 
fm and increase to 38% for Ax = 0.5 fm. At LHC the fluctuations are smaller by 
factor 2.7. Since the average transverse momentum of perturbative partons is of 
GeV order and the longitudinal momentum is even larger, their wavelength is a small 
fraction of 1 fm. The condition (§) is then easily satisfied for the cell size of interest. 
Neglecting of the interaction seems to be reasonable as well, when the partons with 
small 'offshellness' are taken into account. 

Let us summarize our considerations. We have derived the formulas which de- 
scribe how the particle and energy densities fluctuate in a space-time cell of the 
volume Ax 3 At. These formulas, which get a very simple form for massless particles, 
have been applied to estimate the density fluctuations in the quark-gluon plasma. We 
have considered the equilibrium plasma and the anisotropic parton system from the 
early stage of ultrarelativistic heavy-ion collisions at RHIC or LHC. In both cases the 
fluctuations can be large within the reasonable values of the parameters of interest. 
Sizeable density fluctuations are of physical interest when the temporal evolution of 
the plasma system is studied. As shown in 0, the hydrodynamics is then notice- 
ably modified. The density fluctuations seem to be even more important when the 
hadronization is analysed or one considers the processes, such as the J/ip dissociation 
in the plasma [0, which are strongly density dependent. 



I am very grateful to Mark. I. Gorenstein for critical reading of the manuscript. 
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